Memristor for introductory physics 
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Students learn from physics textbooks that there are 
three fundamental two-terminal circuit elements: resis- 
tors, capacitors and inductors. In 2008 May, scientists at 
Hewlett-Packard Laboratories published a paper in Na- 
ture announcing the invention of the fourth element — 
memristor .>i*2 In that paper, Strukov et al presented a 
model to illustrate their invention. Their model provides 
a simple explanation for several puzzling phenomena in 
nanodevices, yet it involves no more mathematics than 
basic algebra and simple integration and differentiation, 
and is comprehensible for college students with knowl- 
edge of first-year physics and calculus. In this note, we 
guide students to derive the analytical solution to the 
equations for the memristor model described in that pa- 
per, and use it to calculate the current responding to a 
sinusoidal voltage. Through this exercise, students will 
realize the nonlinear effects in electric circuit, which is 
ubiquitous in nanoscale electronics. 

The definition of electric current i is the time derivative 
of electric charge q, and according to Faraday's law volt- 
age V is the time derivative of magnetic flux if. In 1971 
Leon Chua noted that there ought to be six mathematical 
relations connecting pairs of the four fundamental circuit 
variables, i, v, q and ipi^ In addition to the definition of 
electric current and Faraday's law, three circuit elements 
connect pairs of the four circuit variables: resistance R 
is the rate of change of voltage with current {R = ^), 
capacitance C is the rate of change of charge with volt- 
age (C = ^), and inductance L is the rate of change of 
magnetic flux with current (L = ^). From symmetry 
arguments, Chua reasoned that there should be a fourth 
fundamental element, which he called a memristor (short 
for memory resistor) M, for a functional relation between 
charge and magnetic flux, M — 
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Resistor 
dv = Rdi 
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Inductor 
dip = L6i 



Capacitor 
dq = Ctiv 



Memristor 
6<p = M6q 



Memristive systems 



FIG. 1: Figure adapted from Nature. 



In the case of linear elements, in which M is a con- 
stant, memristance is identical to resistance. However, 
if M is itself a function of q, yielding a nonlinear circuit 
element, then no combination of resistive, capacitive and 
inductive circuit elements can duplicate the properties 
of a memristor. Almost 40 years after Chua's proposal, 
memristor was found by scientists at HP Labs led by 
R. Stanley Williamsiii^ They suspected that memristive 
phenomenon has been hidden for so long because those 
interested were searching in the wrong places. Although 
the definition of memristor involves magnetic flux, mag- 
netic field does not play an explicit role in the mechanism 
of memristance. As stated in their paper, the mathe- 
matics simply require there to be a nonlinear relationship 
between the integrals of the current and voltage [italics 
added]. The explicit relationship between such integrals 
based on the model by Strukov et all is detailed below. 

After the original proposal of the memristor, Chua and 
Kang generalized the concept to a broader class of sys- 
tems, called memristive systems ,~ defined as 



dx 
'dt 



(1) 



(2) 



Eq. ([T]) looks like Ohm's law, but the generalized resis- 
tance TZ{x) depends upon the internal state x of the de- 
vice. The time derivative of the internal state is a func- 
tion of X and i. Strukov et al were the first to conceive 
a physical model in which x is simply proportional to 
charge q. They designed a device with TZ that can switch 
reversibly between a less conductive OFF state and a 
more conductive ON state, depending on x 



nix) = x{t)Ron + [1 - xit)]Ros 



(3) 



The internal state x{t) is restricted in the interval [0, 1]; 
when x{t) = 0, TZ = Ros, and when x{t) = 1, TZ — Ron- 
The time derivative of x(t) is made to be proportional to 
current 



dx{t) Ro 



dt 



(4) 



in which the parameter (3 has a dimension of magnetic 
fiux (in SI units, magnetic flux is measured in V s, or 
Wb) . The HP scientists produced a semiconductor device 
to satisfy this condition, to be discussed shortly. Eq. ([1]) 
is rewritten as 



v{t) = {x{t)Ron + [1 - x{t)]Ros}i{t) 



(5) 
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Eliminating i{t) from Eq. ([5]) using Eq. (|4]), we can write 
the memristive system as a first-order differential equa- 
tion of x{t). Defining the resistance ratio r — Ron/Rom 
the equation becomes 



v{t) = I3{x{t)+r[l-x{t)]} 



dx{t) 
dt 



(6) 



Using = J vdt, and the property x^ — 5(^2^^), it is 
easy to integrate both sides 



(7) 



where c is a constant of integration determined by the ini- 
tial value of X. This constant contains the past events of 
the input current i and plays an important role in a non- 
linear system. Since x is proportional to charge, Eq. ([7]) 
indicates that the magnetic flux is a quadratic function 
of the charge, and the nonlinear relationship between the 
integrals of the current and voltage for a memristor is 
achieved. 

The memristor invented at HP Labs is made of a thin 
film of titanium dioxide of thickness D sandwiched be- 
tween two platinum contacts. The film has a region 
with a high concentration of dopants having low resis- 
tance Ron, and the remainder has a low dopant concen- 
tration and much higher resistance i?off- The applica- 
tion of an external bias v{t) across the device will move 
the boundary between the two regions by causing the 
charged dopants to drift. Let w{t) be the coordinate of 
the boundary between a less conductive OFF state and a 
more conductive ON state, and /iy the average ion mobil- 
ity (measured in m^ s~^ V^^ in SI miits), Eq. ([4]) written 
in physical parameters is 



1 dw{t) 
D dt 



Ro 



which leads to 



w{t) Ron 



(8) 



(9) 



By comparing Eq. ^ and Eq. ((8|) we identify x = w/D 
and (3 ~ D^/^v- Recall that (3 has a dimension of mag- 
netic flux, thus magnetic fleld is not directly involved in 
this memristor. Eq. ([7]), setting c to zero, becomes 



_ Ronf^V f RoK ^ 



q + Rotiq 



(10) 



For Ros » R on ; 
tained 

M{q) E 



the memristance of this device is ob- 



dip 
dq 



R, 



off 
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fiyRo 



-q 



(11) 



The g-dependent term is due to the nonlinear relationship 
between tp and q. For the flrst time, scientists wrote down 
a formula for the memristance of a device in terms of ma- 
terial and geometrical properties of the device. (An anal- 
ogy is that for parallel plates, capacitance is expressed as 



C = eA/d'^, where e, A and d are permittivity, plate area 
and plate separation, respectively. Similarly resistance 
and inductance of devices can be expressed in terms of 
material and geometrical properties.) Because the coef- 
ficient of q is proportional to memristive behavior 
is increasingly relevant as electronic devices shrink to a 
width of a few nanometers. As semiconductor film of 
thickness D reduces from micrometer scale to nanome- 
ter scale, the nonlinear effect increases by a factor of 
1,000,000. 




FIG. 2: The charge and current of a memristor responding 
to the applied voltage vq sm{ujt) (dashed line in the v-t plot) 
based on Eqs. (I13p and (|14p are shown; the q-tp plot is made 
using Eq. ((Tjl. In these plots, voltage, current, time, magnetic 
flux and charge are expressed in units of vq, vo/Ron, to = 
P/vQ, P and ioto, respectively. The resistance ratio is r = 
Ros/Ro-n = 160; the dimensionless combination of vo/{l3uj) is 
set to be IOO/tt, and (3 (physically / in which D is the 
film thickness and nv the mobility) is chosen to be lO"'^ Wb. 
The i-v plot is produced using the parametric plot procedure 
in graphing software. In the i-v plot, a high-frequency (10a;) 
response is shown, which appears almost as a straight line. 



As soon as researchers discovered the device, it is nat- 
ural to measure the current responding to the applied 
voltage. If the applied voltage in Eq. ^ is 



v{t) = Vq sin(a;t) 



(12) 



then replacing ip = J v{t)dt with — ^ cos(wt) in Eq. (O 
and using the quadratic formula to solve for x{t)^ we have 



xit) 



r2 + 2(r - 1) ^ cos(tjt) c 



1 



(13) 



(We retain only the sign " in "±" so that the solution 
is in the interval [0, 1].) The current i{t) is proportional 
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FIG. 3: This figure is similar to Figured the applied voltage 
is ±vo sin^(tjt) (dashed line in the v-t plot), and resistance ra- 
tio Ros/Ron ~ 380. Notice that while i-v plot in this example 
appears to be complex, the pinched loops feature manifests 
itself vividly. Also the q-ip plot is simple: magnetic flux is a 
single- valued function of charge, see Eq. ([7}. 



We again see the "pinched loops" in the i-v plot. Strukov 
et al reported in their paper that i-v plots shown Fig- 
ures [2] and [3] are regularly observed in their studies of 
titanium dioxide devices. 

The most important lesson for students to learn from 
this note is that the familiar Ohm's law v — iR is merely 
an approximation which is inadequate for a nonlinear 
circuit. Suppose that a scientist, who did not know the 
nature of the titanium dioxide device invented by the HP 
scientists, attempts to characterize this device on the i-v 
plane. If one treats the system as resistive, then based 
on Eqs. (fT^ and (HH), the resistance would be 



R = Rn 



+ 2(r- 1) 



cos{ujt) + c 



(15) 



which is time- varying and frequency-dependent. In 
fact, such an erroneous interpretation of the i-v plot 
is the source of numerous confusions of "current-voltage 
anomalies" reported in unconventional elements jiiS.!^ For 
a memristor, the constitutive relation is a curve in the 
q-(p plane, as depicted in Figures [2] and [3l expressed as 
a single- valued function of the charge q 



if = if{q) 



(16) 



to the time derivative of x{t), which is 



Wo 



r2 + 2(r - 1) 



^ cos(c^i) 



Ron 



(14) 

Our calculation shows that the current is the product 
of v{t) and a function of cos(ujt). Important properties 
of a memristor can be understood from this expression. 
The current i{t) is zero whenever the the applied volt- 
age v{t) is zero, regardless of the state x{t). Because 
the i-v curve must pass through the origin, the i-v Lis- 
sajous figure is double-loop curve — the most salient fea- 
ture of the memristor. At high frequencies, cos(a;<)/ti; 
under the radical sign becomes small, so the nonlinear 
contribution is suppressed. As a result, the i-v hys- 
teresis collapses to a nearly straight line. These fea- 
tures are seen in Figure O The memory effect of a 
memristor is attributed to the constant of integration 
c first appeared in Eq. ([7|). This constant cannot be 
removed by subtracting a constant from x(t), and re- 
mains after differentiating x{t) with time. Therefore, 
memristor response depends on history. As another 
example, with the applied voltage v{t) — iwo sin^ (wt) 
shown in Figure [U the solution is obtained by replacing 
if — ± ^ VQBiv?{ijjt)dt = ±wo[i/2 — cos{ujt) sin(aji)/(2ti))]. 



[This is a rather awkward notation to indicate that ip is 
a function of q.] Differentiating this equation with time 
results in 



dtp dq 
dq dt 



M{q)i 



(17) 



In terms of voltage and current, Eq. assumes the 
form V = M{q)i. Because M{q) varies with the instanta- 
neous value of charge, which book keeps the past events of 
the input current, it has a memory. Without recognizing 
this fact, it is difficult to use v — Ri to explain phenom- 
ena in nonlinear circuits properly. In retrospect, many 
"unusual hysteresis" in i-v plots of various circuits are 
actually due to memristance. With the model of Strukov 
et al, a wide variety of eccentric and complex Lissajous 
figures that have been observed in many thin-film, two- 
terminal devices can now be understood as memristive 
behavior. Such a concept is exceedingly useful in under- 
standing of hysteretic current-voltage behavior observed 
in nanoscale electronic devices, and should be introduced 
to students as early as possible. It is worth mentioning 
that in preparing the paper for Nature, Dr. Williams 
and his colleagues were intentionally writing for curious 
undergraduates.® Students arc encouraged to read the 
original paper for further discussions and potential ap- 
plications. 
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